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LONG-TERM  GOALS 

The  goal  is  to  develop  direct  simulation/physics-based  forward  and  inverse  capabilities  for  radiance 
prediction  in  a  dynamic  ocean  environment.  The  simulation-based  model  includes  and  integrates  all  of 
the  relevant  dynamical  processes  in  the  upper  ocean  surface  boundary  layer  into  a  physics-based 
computational  prediction  capability  for  the  time-dependent  radiative  transfer  (RT). 

OBJECTIVES 

To  develop  physics-based  modeling  and  computational  prediction  and  inverse  capability  for  the  time- 
dependent  underwater  radiative  transfer  incorporating  the  dynamical  processes  on  the  ocean  surface 
and  the  upper  ocean  surface  boundary  layer  (SBL): 

•  Develop  efficacious  direct  computational  capabilities  and  analytical  models  to  predict  the 
radiative  transfer  (RT)  process  of  solar  radiation  in  the  coupled  atmosphere-ocean  system; 
Obtain  understanding,  modeling  and  parameterizations  of  dependencies  of  oceanic  radiance  on 
the  surface  wave  environment. 

•  Develop  direct  simulation  approaches  to  predict  large-scale  realistic  nonlinear  ocean  surface 
wave  fields  for  RT  modeling. 

•  Obtain  quantitative  validation  and  cross-calibrations  with  experimental  field  data. 

•  Investigate  effects  of  dynamic  ocean  surface  waves  on  in-water  solar  radiation  patterns  and 
statistics. 

•  Establish  a  preliminary  framework  for  inverse  sensing  and  reconstruction  of  ocean  surface 
conditions  based  on  underwater  light  field  measurements. 

APPROACH 

We  apply  forward  and  backward  Monte  Carlo  methods  to  simulate  the  radiative  transfer  in  the  coupled 
atmosphere-ocean  system  to  investigate  the  spatial  and  temporal  characteristics  of  the  in-water 
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radiance  and  irradiance  fields.  Both  scalar  and  vector  radiaitve  transfer  simulations  are  developed  to 
consider  the  unpolarized  and  polarized  light,  respectively.  All  our  programs  are  parallelized  and 
optimized  to  run  on  modern  high  performance  computing  (HPC)  platfonns  utilizing  up  to  0(10  ) 
processors.  To  obtain  more  realistic  dynamic  ocean  surface  conditions,  we  extend  and  apply  an  high- 
order  spectral  method  and  LES  capabilities  for  nonlinear  evolution  of  capillary  and  gravity  waves  as 
well  as  their  interactions  with  wind  and  ocean  turbulence  (Xu  et  al  2012).  A  hybrid  method  combining 
Monte  Carlo  and  ray-tracing  is  also  developed  to  characterize  the  radiaitive  transfer  process  in 
inhomogeneous  medium,  such  as  in  strong  ocean  turbulence,  particularly  for  those  with  a  strong 
refractive  index  variation.  More  direct  understanding  of  the  short  term  in-water  light  variability  is 
acquired  by  our  new  physics-based  analytical  model  (we  call  GP)  which  is  able  to  describe  the 
probability  density  function  (PDF)  of  the  downwelling  irradiance  at  the  upper  ocean  (Shen  et  al  2011). 
This  model  takes  into  consideration  the  key  factors  affecting  the  statistics  of  downwelling  irradiance, 
such  as  the  ocean  surface  wave  spectrum  and  the  water  IOPs.  This  model  also  provides  an  asymptotic 
expression  for  the  probability  of  occurrence  of  extreme  values. 

WORK  COMPLETED 

•  Direct  quantitative  comparisons  between  ONR  RaDyO  field  measurements  and  model 
prediction.  We  performed  quantitative  comparisons  by  incorporating  all  the  key  RaDyO 
wave,  IOPs  and  underwater  radiance  measurements  with  combined  wavefield  (reconstruction 
and  prediction)  and  RT  (full  unpolarized  and  polarized  light  field)  modeling  and  closed-form 
analytical  formula. 

•  Data  assimilation  and  preliminary  inverse  modeling.  We  developed  a  framework  and 
obtained  preliminary  results  for  the  inverse  problem  of  deducing  parameters  associated  with  the 
dynamic  ocean  surface  based  on  underwater  light  field  characteristics.  The  immediate  goal  of 
the  inversion  problem  is  to  use  passively  measured  light  field  quantities  including  irradiance, 
radiance  and  polarization  data  in  deterministic  or  statistic  fonns,  to  deduce  ocean  wavefield 
properties  such  as  surface  wind  speed  and  direction,  MSS,  and  wave  frequency  and 
directionality.  Such  inversion  might  be  achieved  utilizing  explicit  formulas  derived  from  the 
analytical  solutions  of  irradiance  and  radiance  statistics  or  spectrum. 

RESULTS 

We  developed  a  3D  Monte  Carlo  (MC)  RT  simulation  capability  for  unpolarized  and  polarized  light 
for  the  atmosphere-ocean  system.  The  model  is  systematically  validated  by  direct  comparisons  with 
existing  theories  and  numerical  model  predictions  and  with  field  data  including  ONR  RaDyo 
measurements.  The  developed  3D  MC  RT  model  was  applied  to  investigate  the  characteristics  of 
polarization  distribution  and  underwater  irradiance  in  various  ocean  surface  environments. 

(1)  Validation  of  models  by  quantitative  comparison  with  field  data:  We  obtained  model-to- 
model  validations  by  comparing  our  simulation  with  the  classic  discrete-ordinate  method  and 
invariant  imbedding  method.  We  made  extensive  cross-calibrations  of  our  simulation 
predictions  with  the  field  data  obtained  in  RaDyO  Santa  Barbara  Channel  (SBC)  and  Hawaii 
experiments  (see,  for  example,  figure  la  and  lb,  for  polarized  light  data  from  Voss  and  IOPs 
data  from  Twardowski).  Validations  of  the  GP  closed-form  model  for  probability  density 
function  of  downelling  irradiance  were  made  by  comparing  the  experiments  data  obtained  by 
Stramski  et  al  and  Monte  Carlo  simulations  (figure  lc). 
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(2)  Effects  of  ocean  surface  roughness  and  nonlinearity  on  underwater  light  fields:  We 
characterized  the  dynamic  ocean  surface  wavefield  in  terms  of  its  roughness  and  nonlinearity. 
The  roughness,  which  is  usually  characterized  by  the  surface  wind  speed  Uio  or  the  mean 
square  slope  (MSS),  has  been  long  recognized  as  the  key  parameter  affecting  the  underwater 
light  field  patterns  and  variability.  However,  quantitative  and  conclusive  studies  of  its  effects 
are  not  clear  yet.  For  both  unpolarized  and  polarized  underwater  light  fields,  we  analyzed  the 
effects  of  the  random  ocean  surface  waves  in  terms  of  the  mean  value  and  higher-order 
moments  of  time-series  light  measurements.  Our  modeling  studies  showed  the  pattern  of  the 
downwelling  irradiance  is  not  affected  by  the  surface  roughness  but  the  variability  decreases 
with  the  wind  speed  and  the  depth  (see  figure  2a).  However,  the  pattern  of  the  degree  of 
polarization  P  decreases  with  wind  speed  (or  MSS)  and  the  depth  (figure  2b);  while  unlike 
irradiance,  the  variability  of  P  increases  with  wind  speed  and  decreases  with  the  depth. 

(3)  Effects  of  wave-induced  light  fluctuations  on  phytoplankton  growth:  Knowledge  of  the  upper 
ocean  light  variability  provides  us  a  new  possibility  to  model  and  explain  the  effects  of  light 
fluctuations  on  the  surface  phytoplankton  growth  characteristics  which  has  been  observed  in 
the  field.  Using  the  analytical  GP  model  we  obtained,  we  showed  that  for  surface 
phytoplankton  exhibiting  photo-inhibition  features  it  is  possible  that  under  strong  light 
illuminations  such  as  in  the  summer  their  short-term  (within  one  day)  growth  rate  can  decrease 
as  the  wind  speed  increases  (figure  3  a).  An  approximate  linear  dependence  of  growth  rate  and 
light  variability  at  the  upper  ocean  were  also  found  (see  figure  3b). 

IMPACT/APPLICATIONS 

The  capability  of  accurate  prediction  of  the  irradiance  transfer  across  ocean  surface  and  in  the  water 
may  enable  the  development  of  a  novel  approach  for  accurate  measurements  of  complex  ocean 
boundary  layer  processes  and  reliable  detection  of  structures/objects  on  or  above  ocean  surface  based 
on  sensed  underwater  irradiance  data. 
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Figure  1.  Validations  of  theoretical  models  by  direct  comparisons  with  experiments:  (a)  Comparisons  of  measured  e- 
veetor  of  sky  polarized  light  fields  (upper)  with  simulated  result  (lower);  (b)  Comparison  of  MC  simulated  in-water 
degree  of  linear  polarization  (DOLP)  and  measured  in-water  DOLP  in  Hawaii  in  the  principal  plane  for  different 
Sun’s  elevations;  (c)  Comparisons  of  PDF  (upper)  and  variance  of  normalized  downwelling  irradiance  y=  Ed/<Ed> 

between  experiments,  MC  simulation  and  analytical  GP  model. 


Figure  2.  Dependence  of  unpolarized/polarized  light  pattern  and  variability  on  ocean  surface  roughness:  (a) 
Variance  of  the  normalized  downwelling  irradiance  x=Ed/<Ed>  as  a  function  of  wind  speed  U10  for  different  depth 

z=-0.86m  ( — ),  -1.7m  ( - ),  -2.85m  (-  •  -)  and  -4.72m  (•••)•  (b)  Dependence  of  degree  of  polarization  ratio  <P>ma,/Psky 

on  the  mean  square  slope  (MSS)  of  the  ocean  surface  at  different  depth  (z=-0.3,  -2.2,  -4.9  and  -7.5  m).  (c)  Variance 
of  the  normalized  degree  of  polarization  P/<P>  at  the  angular  position  of  local  maximum  degree  of  polarization  on 
the  principal  plane  within  the  Snell’s  window  as  a  function  of  wind  speed  Ui0  for  four  different  depths. 


Figure  3.  Effects  of  ocean  wave  roughness  on  surface  phytoplankton  growth  rate:  (a)  Normalized  daily  growth  rate  y 
versus  wind  speed  ulO  for  Anabaena  for  summer  at  different  optical  depths,  (b)  Dependence  ofy  on  scalar 
irradiance  variability  cl  for  Gymnodinium  galatheanum  (green  symbols)  and  Anabaena  (magenta  symbols)  for 
different  optical  depths  Kdz=0.4  (+),  0.6  (A),  0.8  (◄),  1.0  (►),  1.2  (O),  1.4  (*),  1.6  (0)  and  1.8  (□). 
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